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SUMMARY 


The  purpose  of  this  report  was  the  Investigation  of  the  symmetric  modes 
of  motion  of  an  airplane  and  of  the  concepts  which  have  "been  used  to  obtain 
relatively  simpler  expressions  for  the  period  of  these  motions.  The  utility  of 
the  Reeves  Electronic  Analog  Computer  in  this  study  was  investigated. 

Specifically,  the  concepts  concerning  the  symmetric  motions  were  applied 
to  the  general  body  axes  equations  of  motion  in  different  ways.  After  the 
application  of  these  simplifying  assiimptions,  the  modified  equations  were 
solved  for  the  period  of  the  motion  concerned.  Numerical  data  for  a  repre- 
sentative jet  fighter  type  aircraft  at  40,000  feet,  Mach  number  .75^  and  at 
sea  level,  Mach  number  .25  "was  utilized  in  the  solution  by  the  different 
methods.  These  values  for  the  period  were  compared  to  the  true  values  as 
determined  from  the  complete  equations  of  motion. 

The  possibility  of  increasing  the  accuracy  of  the  period  obtained  in 
one  of  the  methods  was  examined. 

A  major  objective  of  this  report  was  the  investigation  of  the  utility 
of  the  analog  computer  in  solving  problems  of  this  type .  Mathematical  treat- 
ment, while  not  difficult,  is  cumbersome  and  slow.  The  analog  computer,  on 
the  other  hand,  after  having  been  supplied  with  information,  achieves  a  solu- 
tion almost  instantaneously.   It  was  desired  to  examine  the  effects  of  changed 
stability  derivatives  and  coupling  between  the  equations  utilizing  the  computer, 
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The  conclusions  reached  as  a  result  of  this  investigation  are: 

1.  The  concepts  concerning  the  short  period  symmetric  motion  of  an  air- 
plane when  applied  to  the  general  equations  provide  an  acceptable 
approximation  to  the  period  of  that  motion. 

2.  Concepts  concerning  the  period  of  the  phugoid  mode,  when  applied, 
provide  accuracies  ranging  from  30  percent  low  to  12  percent  high, 
depending  upon  the  method  of  application.   It  was  concluded,  there- 
fore, that  use  of  the  existing  simplifications  does  not  accurately 
describe  the  period  of  the  true  motion  "because  essential  couplings 
are  neglected. 

3.  The  Reeves  Electronic  Analog  Computer  provided  rapid  accurate  solu- 
tions to  many  sets  of  differential  equations  which  had  to  be  solved. 
The  utility  of  the  computer  in  Investigations  such  as  this  in  unques- 
tionable .  Variation  of  stability  derivatives  or  coupling  between 
equations  could  be  accomplished  by  changing  one  or  several  resistances 
and  the  new  solution  was  immediately  obtainable. 

This  study  was  conducted  at  the  University  of  Michigan,  Ann  Arbor,  Michi- 
gan, during  the  period  February  to  June  1955  under  the  supervision  of  Professors 
M.  A.  Brull  and  J.  D.  Schetzer  of  the  Aeronautical  Engineering  Department. 


INTRODUCTION 

The  advent  of  flexible  high  speed  aircraft  has  somewhat  complicated  the 
problems  which  must  be  solved  by  the  aerodynamic ist  and  the  structures  spec- 
ialist. With  these  complications  has  come  an  increasing  tendency  to  write 
the  equations  of  motion  of  the  airplane  in  the  body  axes  system  of  coordi- 
nates . 

The  purpose  of  this  report  was  the  investigation  of  the  symmetric  modes 
of  motion  of  a  representative  modern  airplane  in  the  body  axes  system  and  in 
particular,  the  investigation  of  the  validity  of  the  simplifying  assumptions 
which  have  been  used  to  obtain  rapid  approximate  solutions  for  the  period  of 
these  motions  o 

The  assumptions  which  have  been  made  are  based  upon  the  observations  of 
Lanchester  who  was  the  first  to  call  the  long  period  motion  the  phugoid. 
Lanchester  stated  that  phugoid  motion  takes  place  at  essentially  constant 
angle  of  attack  and  further  that  the  pitching  motion  is  so  slow  that  the 
aerodynamic  moment  may  be  considered  zero,  or  nearly  so.  When  these  assump- 
tions are  applied  to  the  wind  axes  equations  of  motion  and  additional  terms 
which  do  not  materially  affect  the  period  are  ignored,  an  expression  is 
obtained  for  the  period  of  the  motion  of  .I38  V.  The  same  results  have  been 
obtained  from  considerations  of  kinetic  and  potential  energy  interchange  j 

The  current  literature  does  not  contain  examples  of  the  application  of 
the  phugoid  simplifying  assumptions  in  the  body  axes  equations  of  motion. 


The  short  period  mode  of  motion  is  defined  as  a  motion  in  pitch  of  very 
short  duration  which  takes  place  as  the  angle  of  attack, which  has  heen  changed 
due  to  a  disturbance,  returns  to  the  equilibrium  value.  The  motion  is  so  rapid 
that  it  is  assumed  that  the  forward  speed  does  not  change  and  consequently, 
the  motion  is  independent  of  changes  in  perturbations  in  the  forward  speed. 


SYMBOLS 

Vp  aircraft  velocity  ft/sec 

X  aerodynamic  force  in  X  "body  axis  direction  lbs 

Z  aerodynamic  force  in  Z  body  axis  direction  lbs 

M  aerodynamic  moment  about  the  Y  body  axis  ft  lbs 

u  perturbation  velocity  in  X  body  axis  direction  ft/sec 

w  pertiirbation  velocity  in  Z  body  axis  direction  ft/sec 

q  perturbation  velocity  in  Y  body  axis  rad/sec 

dot  over  a  symbol  -  derivative  of  symbol  with  respect  to  time 

X  operator  for  derivative  with  respect  to  time,  — 

X•^   dimensional  derivative  of  the  aerodynamic  force  in  the  X  body  axis 

^X 
direction  with  respect  to  the  perturbation  velocity  u;  ^  lb  sec/ft 

Xw   1^  lb  sec/ft 
ow 

Z^^   l|  1^  sec/ft 

Z^   ^  lb  sec/ft 
^w 

Zn   M  lb  sec/ft 

M^   ^  f t  lb  sec/ft 
du 

M^   ^  f t  lb  sec/ft 
^w 

M^   M  f t  lb  sec/ft 

Mq   ^  ft  lb  sec/radian  ' 
dq 


QL  angle  of  attack  of  M.A.C.  relative  to  wind  -  degrees 

6i  angle  of  pitch  of  X  axis  measured  from  the  horizon  -  degrees 

Ui  steady  state  lihear  velocity  along  X  axis  -  ft/sec 

Wi  steady  state  linear  velocity  along  Z  axis  -  ft/sec 

m  mass  in  slugs 

lyy  moment  of  inertia  about  Y  body  axis  -  slug  ft 

M  Mach  number 

Cl  Lift  coefficient 

F^  Sum  of  forces  in  X  body  axis  direction 


SYMMETRIC  EQUATIONS  CfF  MOTION 
(BODY  -  AXES  SYSTEM) 

The  six  equations  of  motion  "break  down  into  two  sets  of  three  equations. 
One  set  describes  asymmetric  motion.   The  other  set  describes  symmetric  motion, 
The  linearized  equations  for  symmetric  motion^  neglecting  X^,  in  "body  axes 
coordinates  are: 

F^  =  0  =  [X^^-mXju  +  [X^jw  +  [-mWi-2i-£°£-®i]  q 
F^  =  0  =  [Z^]u  +  [Z^-mXjw  +  [Zq-fmUi-^  ^^^  ^3  ]q 
My  =  0  =  [M^Ju  +  [Mv^l^^Ow  +     [M^  -  lyy  X.  ]  q 

The  determinant  of  the  coefficients  is  set  equal  to  zero  and  the  stabi- 
lity quart ic  is  obtained 

AX*  +  BX^  +  CX^  +  DX  +  E  =  0 
where 


A  = 

1 

B  = 

m 

^1         ^a              Ma 

C   = 

"^■^yy                       -'-yy 

g  sin  9i 


Zq+mUi  1 

+  ^Fl^T"  ( Vu-^)  +  ^  (XuZ^-^uXw) 

^               s  sin  ©1  Wt 

D  =  -^^   (Z^^Xv-X^Zw)  -  ^ ^(M^Zu-M^,Zv,)  +  -^   (M^Z^,-M^Z^) 

(Zq+mUi)  X^M^  g  cos  01  ,,,  ^   ,,  , 

miyy  mlyy 


The  equations  of  motion  are  developed  in  Ref .  1. 


When  the  analog  computer  is  to  be  utilized,  it  is  convenient  to  rewrite 
the  equations  of  motion  as  follows: 

ml  =  XyU  +  X^w  -  m  W^q-Mg  cos  9i  /q 
ra^  =  Zyu  +  Z^w  +  [Zq+mUi]q  -  mg  sin  @xjq 
lyyq  =  I^u  +  M^^  +  1^^  +  Mqq. 

The  computer  circuit  for  the  solution  of  these  equations  is  included  as 
Fig.  1.  The  airplane  niomerical  data  supplied  to  the  circuit  is  controlled 
hy  positioning  the  potentiometers.  Each  potentiometer  controls  a  specific 
airplane  parameter.  These  are: 

1.  m  6.  X^  11-  ^ 

2.  X^  7.  mWi 
5.  m  8.  mg  sin  9i 

5.    lyy  10.   M^ 


12. 

\ 

13. 

Zq  +  mUi 

Ik. 

Mq 

15. 

mg  cos  @2. 

APPLICATION  OF  CLASSICAL  ASSUMPTIONS 
FOR  SIMPLIFICATION  OF  THE  SOLUTIONS 
TO  THE  MODES  OF  MOTION 

The  roots  of  the  stahility  quartic  yield  the  characteristics  of  the 
short  period  and  phugoid  motions.  Various  simplifying  assumptions  con- 
cerning these  motions  have  "been  introduced  in  the  past.  Accuracy  obtain- 
able after  the  assumptions  have  been  applied  varies.   Normally,  the  short 
period  motion  is  so  heavily  damped  and  so  rapid  that  it  is  of  little  con- 
sequence .   It  does  become  important  if  the  period  is  within  the  normal 
human  response  lag  and  not  heavily  damped,  primarily  because  a  pilot  may 
inadvertently  reinforce  the  motion  when  he  attempts  to  correct  it .  The 
short  period  mode  is  important  when  platform  stability  is  a  consideration. 
The  phugoid  motion  is  a  long  period  slow  oscillation  with  weak  damping  and 
of  such  little  consequence  that  it  may  be  negatively  damped  and  not  render 
an  airplane  unacceptable.  Because  of  the  apparent  unimportance  of  these 
modes  in  piloted  subsonic  aircraft,  approximations  to  their  period  and 
damping  are  justified  and  desirable. 

However,  the  importance  of  these  modes  of  motion  in  pilotless  missiles 
at  very  high  speeds  has  not  been  completely  investigated.  When  a  relatively 
long  slender  missile  encounters  severe  turbulence  or  a  sharp  gust,  the  pos- 
sibilities seem  to  be  present  for  considerable  control  difficulties.  These 
motions  could  intensify  the  difficulties . 

It  follows  then  that  it  is  important  to  have  a  method  for  the  rapid, 
accurate  determination  of  the  symmetric  modes . 


SHORT  PERIOD  METHOD 

The  short  period  motion  is  obtained  by  assiiming  that  the  motion  proceeds 
at  constant  forward  velocity  (u  =  0).  This  follows  if  it  is  reasoned  that  by 
the  time  the  forward  velocity  has  changed,  the  motion  has  damped  out  and  has 
ceased  to  exist .  The  assumption  permits  a  reduction  of  the  equations  of 
motion  to  essentially  a  coupling  between  the  perturbations  w  and  q.  The 
equations  of  motion  become: 


F^  =  0  =  [Z^-m\]w   +   [Z^+mUi-  Si-^iL©!], 


My  =  0  =  [M^+M^\Jw 


+      [Mq  -  lyy  ^.]q 


Solution  of  these  equations  as  written  involves  the  solution  of  a  cubic  . 
The  equation  becomes  a  quadratic  if  the  perturbation  gravity  force  is  assiimed 
very  small  and  ignored.  The  quadratic  is: 

AX^  +  B\  +  C  =  0 
where : 


A  =  1 


B  =  - 


lyy     m       miyy 


;Zq+mUi) 


miyy  ^lyy  ^ 


The  utility  of  the  analog  computer  is  demonstrated  by  the  fact  that  the 
computer  circuit  for  the  short  period  method  is  achieved  by  uncoupling  poten- 
tiometers 11  and  12  in  Fig.  1.  This  uncouples  and  removes  from  the  circuit 
all  contribution  to  the  short  period  of  the  u  perturbation  velocity  and  the 
solution  is  obtained  using  the  reduced  Z  force  and  moment  equations . 

Ref .  2  shows  a  simplified  short  period  method  in  the  wind  axes  system. 
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PHUGOID  METHOD 

The  long  period  phugoid  motion  is  usually  obtained  by  applying  Lan- 
chester's  original  assumptions.  Ref .  3  states  that  in  a  phugoid  motion  the 
angular  pitching  movement  of  the  aircraft  is  at  all  times  very  slow  and  this 
implies  that  the  aerodynamic  pitching  moment  is  nearly  zero.  Hence,  the 
angle  of  attack  is  at  all  times  very  nearly  that  corresponding  to  zero 
pitching  moment.  This  implies  nearly  constant  angle  of  attack.  Ref.  5 
further  assiimes  that  o;  and  Cj^  is  strictly  constant. 

This  investigation  revealed  that  a  choice  is  available  for  the  method 
of  enforcing  zero  change  in  angle  of  attack  in  the  body  axes  equations. 

By  definition: 

OL  =   arctan  — ^ 
Ui 

if  dWi  =  perturbation  velocity  in  Z  direction  =  w 
and  dUi-=  perturbation  velocity  in  X  direction  =  u 
then 

Aq;  =  -—  w  +  -—  u 

.  ,   w  cos  q;  -  u  sin  OL 
ha  -  


PHUGOID  METHOD  I 

The  expression  for  the  perturbation  angle  of  attack  is  usually 
simplified  by  assuming  the  equilibrium  angle  of  attack  is  small,  and 


because  of  this  assumption,  that 

Aa  ^  J;L 
Vp 

By  this  method,  the  perturbation  a   is  related  directly  to  the  pertur- 
bation velocity  w. 

When  Lanchester's  assumptions  are  applied  to  the  equations  of  motion 
based  on  this  near  equality  for  La,   and  Ace  =  w  =  0,  the  stability  quart ic 
reduces  to  a  quadratic  which  gives  an  approximate  solution.   The  quadratic 
is  obtained  by  ignoring  the  moment  equation  and  deleting  all  terms  in  w  in 
the  force  equations.  Two  equations  remain  in  the  variables  u  and  q.  These 


ares 


Fx  =  0  =  [X^-m>,]u   -    [mWi+2S_£O£_01]q 
F^  =  0  =    [Z^]  u  +  [Zq+mUi-M_2i2_®iJq 

The  determinant  of  the  coefficients  set  equal  to  zero  yields: 

AX^  +  B\  +  C  =  0 
Where: 

A  =  [Zq+mUi] 

Z. 
B  =  -[X^  (^  +  Uj)  +  mg  sin  61+  Z^,Wi] 

C  =  g  [X^  sin  61-  Z^^.  cos  e^j 

When  the  analog  computer  is  utilized,  the  two  equations  are  rewritten 

as  follows: 

mu  =  X^il  -  mW^q  -  mg  cos  8iq 

(Zq+mU3)q  =  mg  sin  8iq  -  Z  li 

The  computer  circuit  is  shown  in  Fig.  2.   Specific  airplane  parameters 

are  controlled  by  the  potentiometers  as  follows: 

1.   Zq  +  mlJi  2.  mg  sin  9i 
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5.  m  ^-     \ 

?•  X^  8.  mg  cos  9i 

9.  mWi 
Switch  in  Fig.  2  is  open. 

This  method  of  simplification  produces  a  period  which  is  greater  than 
the  true  period. 

PHQGOID  METHOD  II 

A  second  method  for  the  reduction  for  the  complete  equations  to  obtain 
the  period  of  the  phugoid  mode  was  developed  by  the  writers  during  this  in- 
vestigation. This  method  involves  the  application  of  the  identical  simpli- 
fying assumptions  of  method  I,  but  the  complete  equality 

Aa  =  ^   ^Q^  Q!  -  u  sin  PL   _  0 

was  used  to  enforce  no  change  in  angle  of  attack.  Then  w  =  u  tan  CH  was 
substituted  in  the  equations  and  the  third  equation  was  ignored.  The  three 
equations  reduce  to  two  equations  which  are: 

F„  =  0  =  [X,-m\  +  X^  tan  a]u   -  [m  Wi  +-SS_2^_§i]q 

■^  A. 

F^  =  0  =  [Zu+Z^tana-m\tana]u  +  [Zq+mUi-°^  ^^^®l]q 
The  determinant  of  the  coefficients  set  equal  to  zero  yields: 

AX^  +  B\  +  C  =  0 


where : 


A  =  m  [Ui+Wi  tan  a   +  ^] 

m 

B  =  -  [(Zq+mUi)(--iif  :^  tan  a)  +  mg  (sin  01  -  cos  8i  tan  a) 
+  Wi  (VZ^  tan  a)] 
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C  =  g  [sin  e^  (X^  tan  o:  +  X^)  -  cos  8i  {7^7^   tan  a)] 
The  equations  are  rewritten  as  follows  when  the  analog  computer  is  used: 
mx     =   (X^+X^  tan  a)  U  -  m  W^q  -  mg  cos  9i  q 
(Zq+m  Ui)q  =  mg  sin  81  q  -  (Zu+Zv  tan  a)  u  +  mil  tan  05 
The  computer  circuit  is  shown  in  Fig.  2.   Specific  airplane  parameters 
are  controlled  "by  the  potentiometers  as  follows: 

1.   Zq  +  mUi         6.   Z^  +  Zv7  tan  a 


mg  sin  81 

m 

m  tan  OC 

m 


7.  X;;^  +  X^  tan  a 

8.  mg  cos  81 

9.  mWi 


Switch  in  Fig.  2  is  closed. 

This  method  of  simplification  produces  a  period  which  is  less  than  the 
true  period. 

Solution  "by  Method  II  provides  a  period  which  is  exactly  equal  to  .I38 

V 
V-r,,  arrived  at  in  wind  axes  in  Ref.  2,  and  is  exactly  equal  to  k.kk  -^     de- 
^  g 

rived  in  Ref.  3.  Thus  Method  II  is  the  true  application  of  Lanchester's 
Basic  Assiomptions  in  "body  axes.  These  assumptions  are  based  on  the  pure 
interchange  of  kinetic  and  potential  energy. 

Ref.  k   states,  "To  secure  a  lift  proportional  to  the  square  of  the  speed, 
the  wings  must  exert  a  constant  lift  coefficient,  so  that  their  incidence  must 
be  constant.  This  implies  that  the  aeroplane  is  statically  stable  as  regards 
pitching,  that  it  has  a  negligible  moment  of  inertia  about  the  pitching  axis, 
and  the  length  between  wings  and  tail  is  small  compared  to  the  radii  of  cur- 
vature of  the  flight  path." 
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Method  II  does  not  substantiate  this  statement  since  substitution  of 
w  =  u  tan  a,    (for  the  condition  Act  =  0)  in  the  moment  equation  causes  M^^ 
and  M^  to  be  grouped  in  the  same  coefficient  and  these  terms  exactly  cancel 
each  other  because  they  become  equal  in  magnitude  and  opposite  in  sign. 
(  \}\x\   =    \^   "tan  Q;  I )  .  This  leaves  the  moment  equation  without  a  static  sta- 
bility contribution  and  produces  the  same  effect  as  zero  static  stability. 
In  fact,  solutions  of  pairs  of  the  simplified  equations  with  the  moment 
equation  included  as  one  of  the  pair  show  the  motion  to  be  aperiodic .  This 
is  the  motion  to  be  expected  from  an  airplane  with  zero  static  stability. 
The  condition  of  zero  static  stability  may  be  solved  quickly  on  the  analog 
computer  by  positioning  the  two  potentiometers  which  control  M^  and  M^  on 
the  total  airplane  circuit  to  zero.  Computer  solutions  for  this  case  show 
aperiodic  motion. 

It  is  apparent  that  the  Lanchester  idealization  enforces  neglect  of  the 
moment  equation  because  as  far  as  the  phugoid  mode  is  concerned,  the  moment 
equation  becomes  inconsistent  with  the  two  force  equations.  This  requires 
that  the  two  force  equations  be  used  to  obtain  the  solution  for  the  phugoid 
mode  and  that  an  error  in  period  must  be  accepted. 

MODIFIED  PHUGOID  METHOD  II 

It  should  be  noted  that  Aa  =     ^     =     0,  if  applied  in  the  body  axes 

Vp 

equations,  produces  a  period  which  is  larger  than  the  true  period;  whereas. 
Method  II  produces  a  period  which  is  smaller  than  the  true  period. 

Comparison  of  the  determinantal  coefficients  resulting  from  both  Aa 
approximations  was  made  and  the  terms  in  each  of  the  coefficients  which 
largely  effected  the  period  were  examined.   It  was  noted  that  the  gravity 
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contribution  in  the  coefficient  multiplying  \  to  the  first  power  (Method  II) 
was  completely  destroyed  by  tan  CX. 

[mg  (sin  61  -  cos  9i  tan  a)]      a  =  9i 

The  classical  explanation  of  phugoid  motion  is  based  on  the  interchange 
of  kinetic  and  potential  energy.  The  complete  removal  of  all  effects  of 
gravity,  therefore,  seems  illogical.   In  short,  if  the  factor  tan  a  is  re- 
moved from  the  term  mk   tan  o:  (second  equation),  the  gravity  effects  are 
restored  and  the  solution  for  the  period  of  the  phugoid  is  accurate  to  with- 
in one  percent .  The  resultant  period  is  nearly  the  actual  value . 

This  result  is  an  experimental  result  arrived  at  on  physical  grounds 
alone  and  perhaps  may  not  be  applied  with  such  accuracy  in  general  because 
it  cannot  be  mathematically  justified.  As  shown  above,  the  idealizations 
of  Lanchester  as  generally  applied  are  also  not  mathematically  Justified. 
(When  Aa  =  0,  the  airplane  is  not  simultaneously  statically  stable  in  the 
equations  of  motion.)   If  the  additional  reasoning  concerning  the  gravity 
term  is  accepted  as  described  on  physical  grounds,  the  coefficients  of  the 
quadratic  describing  the  phugoid  mode  are: 

A  =  m  [Ui+Wi  +  ^] 
B  =  -  [(Zq+mUi)(^+^  tan  a)  +  mg  (sin  61  -  cos  61) 

+  Wi  (Zru+Zv,  tan  a)] 
C  =  g  [sin  81  (Xu+X^  tan  a)   -   cos  81  (Z,u+2^  "tan  a)] 

The  computer  solution  for  Modified  Method  II  is  obtained  by  charging  the 
setting  of  potentiometer  4  in  Fig.  2. (switch  closed)  to  remove  the  factor  tan  cc. 

Numerical  examples  for  Method  I,  Method  II,  and  Modified  Method  II  are 
included  in  this  report.  Method  I  is  12  percent  high.  Method  II  is  seven 
percent  low  and  Modified  Method  II  is  within  one  percent,  all  referred  to 
the  actual  period. 
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PHUGOID  METHOD  III 

One  other  method  for  an  approximate  solution  for  the  phugoid  mode  has 
"been  used.  This  method  makes  use  of  the  fact  that  phugoid  motion  is  essen- 
tially a  coupling  hetween  u  and  /q.  The  method  of  simplification  used  is  to 
delete  the  w  coupling  terms  in  the  force  and  moment  equations  thus  uncoupling 
the  Z  force  equation  from  the  system.  The  X  force  equation  and  the  moment 
equation  thus  reduced  are  solved  simultaneously.  The  equations  are: 

Fx  =  0  =  [Xu-mXju  -  [mWi+2S_£0L®i]q 

Mj,  =  0  =   [  M^  ]u  +  [  Mq  -  lyy  \   ]q 
The  resulting  cubic  equation  is: 

AK^  +  D\2  +  CX  +D  =  0 


where : 


A  =  1 


B 


m  Mq  +  lyy  X^ 


°^  V 


fXu  Mq  +  m  Wi  Mu"| 

L     m  lyy        J 

L    Iw    J 


C  = 

■-yy 


■■yy 

This  method  produces  an  error  in  the  period  of  25  to  30  percent . 

Solutions  on  the  Reeves  Electronic  Analog  Computer  show  this  error . 
The  computer  solution  is  obtained  by  deleting  the  connections  to  potentio- 
meters 6,   9,   and  10  in  Fig.  1.  This  effectively  uncouples  the  Z  force  equa- 
tion. The  most  powerful  term  in  reducing  the  period 
was  determined  to  be  [M^jw.  This  was  shown  by  step  by  step  replacement  of 
the  deleted  terms.  When  [Mv/]w  alone  was  reconnected,  the  period  was  restored 
to  the  correct  value . 
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One  other  idealization  was  applied  to  the  equations  of  motion  in  an 
attempt  to  find  a  simple  accurate  representation  of  the  phugoid  mode.  Because 
the  phugoid  is  described  as  a  long  slow  oscillation  about  a  datum  line,  it 
seemed  reasonable  to  assume  that  motions  of  this  type  would  be  characterized 
by  very  small  accelerations  in  the  direction  normal  to  the  flight  path.  This 
may  be  expressed  as: 

w  -  Ui  q  =  0 
This  equation  was  solved  for  q  and  substitution  for  q  was  made  in  the  equations 
of  motion.   Pairs  of  the  equation  were  solved  with  varying  results.   Solution 
of  the  two  force  equations  produced  aperiodic  motion;  solution  of  the  Z  force 
equation  and  the  moment  equation  produced  a  short  period  mode  frequency;  and 
solution  of  the  X  force  equation  and  the  moment  equation  produced  a  cubic . 
This  cubic  was  not  solved  but  known  accurate  values  of  the  phugoid  roots  were 
substituted  to  ascertain  whether  or  not  the  correct  phugoid  mode  was  contained 
in  the  cubic  .   It  was  found  that  the  roots  did  not  satisfy  the  equation  even 
remotely  and  the  method  was  discarded. 

It  seemed  logical  to  assiome  further  that  if  the  normal  accelerations  are 
small  or  zero,  that  the  condition  of  zero  angle  of  attack  change  would  occur 
simultaneously.  When  this  additional  condition  was  imposed,  the  equations  of 
motion  reduced  to  three  equations  in  u  but  all  the  solutions  indicated  aperio- 
dic motion. 

NUMERICAL  RESULTS  AND  DISCUSSION 

The  following  examples  are  selected  to  facilitate  comparison  of  the 
results  obtainable  by  the  several  methods  discussed  above.  The  airplane 
upon  which  this  is  based  is  a  representative  jet  fighter  type  aircraft. 
Two  case  are  presented:   first,  the  case  of  Mach  number  0.75  at  40,00  feet; 
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second,  the  case  of  Mach  niimber  0.25  at  sea  level.   Noinml  mathematical 
solutions  are  presented  and  then  are  verified  by  Reeves  Analog  Computer 
solutions  for  the  first  case  „  The  second  case  will  be  presented  as  analog 
computer  solutions  only. 

CASE  I 
AIRPLANE  DATA  AND  VELOCITIES 

S  =  600  sq  ft  Xu  =  -5.19  Ih  sec/ft 

m  =  1180  slugs  X^  =  1^8.8  lb  sec/ft 

lyy  =  80,000  slug  ft2  Z^  =  -70.8  lb  sec/ft 

h  =  i4-0,000  ft  altitude  Z^  =  -807  lb  sec/ft 

a  =  2.i45  degrees  =  81  U^  =   -637  lb  sec 

Vp  =  728  ft/sec  M^  =  10.08  lb  sec2 

Ui  =  728  ft /sec  M^  =  27.5  lb  sec 

Wi  =  31.5  ft/sec  Zq  =-k.9k   x  10^  lb  sec 

M  =  0.75  Mq  =  -128. i+  X  103  It  ft  sec 

Using  this  data,  the  stability  quart ic  coefficients  are: 

A  =  1 

B  =  2.2051+ 

c  =  6.8699 

D  =  .01+604 
E  =  .02291 

By  Bairstow's  approximate  factorization  and  one  iteration  for  the  phugoid 

roots,  the  quadratic  is: 

\^  +  .0056i+\  +  .00332  =  0 

The  period  of  the  phugoid  is  109  seconds  and  the  time  to  damp  to  half 

amplitude  is  21+6  seconds , 
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The  short  period  roots  were  obtained  from  the  remaining  factor  of  the 
quartic: 

X^  +  2.199T6\  +  6.85^^19  =  0 
The  period  of  the  short  period  mode  is  2.65  seconds  and  time  to  damp  to  half 
amplitude  is  O.63  seconds. 

The  complete  solution  was  obtained  by  using  the  Reeves  Electronic  Ana- 
log Computer  and  these  results  are  shown  in  Fig.  3.  The  values  obtained  by 
the  computer  method  and  the  stability  quartic  method  are  seen  to  be  in  excel- 
lent agreement  and  these  values  will  be  used  as  the  basis  for  comparison  of 
the  various  methods  employed  for  simplification. 

SHORT  PERIOD  METHOD 

It  is  assumed  that  the  short  period  mode  is  essentially  a  coupling  be- 
tween the  perturbations  w  and  q  alone.  By  this  method,  the  equation  is: 

A\2  +  B\  +  C  =  0 
where : 

A  =  1 

B  =  2.201 

C  =  6.8i^ 
From  this,  the  period  is  2.66  seconds  and  the  time  to  damp  to  half  amplitude 
is  0.63  seconds.  This  is  in  excellent  agreement  with  the  results  obtained 
from  the  quartic.  The  computer  solution  is  shown  in  Fig.  k. 

PHUGOID  METHOD  I 

The  second  method  to  be  compared  is  the  previously  discussed  phugoid 

simplification  in  which  Ace  was  considered  to  be  equal  to  JL  .  The  quadratic 

Vp 
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equation  for  this  method  is: 

A?\.2  +  B\  +  C  =  0 
where : 

A  =  8.5306 

B  =  .0^356 

c  =  .02275 

The  period  hy  this  method  is  122  seconds  and  the  time  to  one-half  amplitude 
is  273  seconds.  The  period  is  seen  to  be  about  12  percent  in  error.  The 
computer  solution  is  shown  in  Fig.  5. 

PHUGOID  METHOD  II 

The  third  method  to  be  compared  is  the  method  which  requires  that  when 
Aa  =  0,  w  =  u  tan  a.  The  quadratic  equation  which  results  is: 

AX^  +  BK  +  C     =     0 
where : 

A  =  8.56 
B  =  .0555 

c  =  .031+ 

The  method  gives  a  period  of  101  seconds  and  the  time  to  damp  to  half  ampli- 
tude is  213  seconds.  The  computer  solution  is  shown  in  Fig.  6  for  further 
comparison.   It  is  seen  that  by  this  method  the  period  is  about  eight  percent 
too  low. 

MODIFIED  PHUGOID  METHOD  II 

The  fourth  method  to  be  compared  is  the  method  evolved  in  this  report. 
This  method  requires  that  when  Aa  =  0,  w  =  u  tan  ex.     The  quadratic  equation 
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is: 

AX2  +  B\  +  C  =  0 
where : 

A  =  8„92 

B  =  MQ 

c  =  .03U 

The  period  by  this  method  is  110  seconds  and  the  time  to  damp  to  half  ampli- 
tude is  29.6  seconds.  The  computer  solution  is  shown  in  Fig.  7  for  fiirther 
comparison.   It  is  seen  that  by  this  method  the  period  is  accurate  to  less 
than  one  percent . 

PHUGOID  METHOD  III 

The  fifth  and  final  numerical  solution  is  the  last  phugoid  case  where 
w  is  considered  to  be  zero  and  the  second  force  equation  is  ignored.  This 
method  produces  the  cubic  equation: 

A\^  +  BX^   +  C\  +  D  =  0 
where : 

A  =  1 

B  =  1.60939 

c  =  .01788 
c  =  .01106 

This  produces  one  negative  real  root  and  a  complex  pair.  The  period  by  this 
method'  is  75.7  seconds  and  the  time  to  damp  to  half  amplitude  is  203  seconds 
The  computer  solution  is  shown  in  Fig.  8  for  further  comparison.  The  error 
by  this  method  is  about  30  percent . 
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CASE  II 
AIRPLANE  DATA  AND  VELOCITIES 

S  =  600  sq  ft  Xu  =  -6.3  lb  sec/ft 

m  =  1180  slugs  X^  =  165  lb  sec/ft 

iyy  =  80,000  slug  ft2  Zu  =  -162.6  lb  sec/ft 

h  =  sea  level  Z^  =  -98O  lb  sec/ft 

a  =  6.5  degrees  =  @i  ^  =  -13^5  lb  sec 

Vp  =  279  ft/sec  M^  =  -102.7  lb  sec^ 

Ui  =  279  ft/sec  M^  =  153.7  lb  sec 

Wi  =  31.9  ft/sec  Zq  =  -^325  lb  sec 

M  =  .25  Mq  =  -115000  lb  ft  sec 

Computer  solutions  for  the  second  case,  Mach  0.25  at  sea  level  are 

shown  in  Figs .  9  throu&h  I3  and  the  results  are  seen  to  parallel  Case  I . 

CONCLUSIONS 

Comparison  of  the  numerical  results  and  computer  solutions  of  the 
several  methods  presented  lead  to  the  following  conclusions: 

1.  The  concepts  concerning  the  short  period  symmetric  motion  of  an 
airplane  when  applied  to  the  general  equations  provide  an  acceptable 
approximation  to  the  period  of  that  motion. 

2.  Concepts  concerning  the  period  of  the  phugoid  mode,  when  applied, 
provide  accuracies  ranging  from  30  percent  low  to  12  percent  high, 
depending  upon  the  method  of  application.   It  was  concluded,  there- 
fore, that  use  of  the  existing  simplifications  does  not  accurately 
describe  the  period  of  the  true  motion  because  essential  couplings 
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are  neglected , 
3.  The  Reeves  Electronic  Analog  Computer  provided  rapid  accurate  solu- 
tions to  many  sets  of  differential  equations  which  had  to  be  solved. 
The  utility  of  the  computer  in  investigations  such  as  this  is  unques- 
tionable. Variation  of  stability  derivatives  or  coupling  between 
equations  could  be  accomplished  by  changing  one  or  several  resistances 
and  the  new  solution  was  immediately  obtainable . 
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